Correlated diffusion of colloidal particles near a liquid-liquid interface 
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Abstract - Optical microscopy and multi-particle tracking are used to investigate the cross- 
correlated diffusion of quasi two-dimensional (2D) colloidal particles near an oil-water interface. 
It is shown that the effect of the interface on correlated diffusion is asymmetric. Along the line 
joining the centers of particles, the amplitude of correlated diffusion coefficient Dy (r) is enhanced 
by the interface, while the decay rate of Dy (r) is hardly affected. At the direction perpendicular to 
the line, the decay rate of D^{r) is enhanced at short inter-particle separation r. This enhancing 
effect fades at the long r. In addition, both D\\ (r) and D±{r) are independent of the colloidal area 
fraction n at long r, which indicates that the hydrodynamic interactions (His) among the particles 
are dominated by that through the surrounding fluid at this region. However, at short r, D±{r) is 
dependent on n, which suggests the His are more contributed from the 2D particle monolayer self. 
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I. INTRODUCTION 



Much attention has been attracted to the dynamic be- 
havior of confined colloidal suspensions recently [l|-|6| . In 
real circumstances, particles are usually spatially con- 
fined in special environments, such as microfluidic de- 
vices, porous media, fluid interface and cell membrane 
[THllj. The dynamical behaviors of colloids in spatially 
confined environments are more complicated than that 
in unbounded three-dimensional (3D) fluid bulk. In the 
unbounded 3D bulk suspensions, the longitudinal and 
transverse correlated diffusion coefficient D\\ and D± are 
well known as D|| a l/r, D± cc l/r and Z?|| = 2D± 
[13 |. For the particles confined near the solid wall, the 
strength of hydrodynamic interactions (His) decay with 
inter-particle separation r as 1/r^ [J, [ij, ^^. In addi- 
tion, spatial symmetry break due to the boundary con- 
ditions brings asymmetry feature to cross correlated dif- 
fusion of particles 0, I12h15| . For the particles on the 
air-water interface or a viscous membrane, both experi- 
mental [1^ and theoretical [1^ studies show Z)|| oc l/r 
and D± oc 1/r^. Oppenheimer's et al. calculation |17| 
showed that at the membrane adjacent to a solid wall, 
the correlated diffusion of particles is a function of the 
distance between the membrane and solid wall. 



The influence of different boundary conditions could be 
directly reflected on diffusion behavior of particles [isl - 
|23 |. Many studies base on solid wall condition, a non-slip 
boundary, which could cut off the fluid fleld in its vicin- 
ity. For a fluid-fluid interface (say a soft wall), there is 
a slip boundary which partially transforms the flow field 
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around. Compared with extensive investigations on the 
solid wall's effects, there are few experimental studies 
about the infiuences of the fiuid-fluid interface on parti- 
cles' correlated diffusion. 

In this paper, we report the experimental measurement 
of the cross-correlated diffusion of colloidal particles near 
the water-decahydronaphthalene (decalin) interface. The 
particle monolayer is similar to a viscous membrane. The 
correlated diffusion of particle pair shows asymmetric be- 
havior in two normal directions. Along the line joining 
the centers of particles, the decay rate of correlated dif- 
fusion is constant. At the perpendicular direction, the 
decay rate increases with the particle separation. In both 
directions, the correlated diffusion is independent of the 
area fraction of particles for long particles separation. 
The influence of fluid-fluid interface on the dynamic be- 
havior of the colloidal monolayer tends to be saturated, 
when the particle separation is much larger than the dis- 
tance between the interface and monolayer. These phe- 
nomena reveal weights of His through different paths. 

The remainder of the paper is organized as follows. 
Section II is devoted to the description of the experi- 
mental setup. Experimental results and discussions are 
presented in Sec. Ill . Finally, the work is summarized 
in Sec. IV. 



II. EXPERIMENTAL SETUP 

Two kinds of colloidal particle with the same diameter 
d —2.0vfj,m are used in the experiment. The flrst one is 
polystyrene (PS) latex spheres purchased from Invitro- 
gen with sulfate and high density of carboxyl functional 
groups on the surface. The other one is silica spheres 
purchased from Bangs, which have anionic SiO~ groups 
on the surface when dispersed in water. The density of 



the silica spheres are about 1.9 times larger than that of 
PS. Both the polystyrene and silica spheres are the repre- 
sentative of charged particles commonly used in colloidal 
science. 

The experimental setup is similar to that described in 
Ref. 23. The sample cell is made from a stainless steel 
disk with an inner container in which there is a hole with 
the diameter 8.3 mm. The bottom of the hole is sealed 
with a 0.1 m,m thick glass cover slip, which also serves as 
an optical window. We first fill the hole to the top edge 
with deionized water of 18.2 Mil ■ cm, in which the par- 
ticles cleaned 7 times by centrifugation were suspended. 
And then we add decalin (a mixture of cis and trans with 
density 0.896 g/cm.^) to the top of the water, filling the 
entire height of the inner container. The solvent decalin 
was purchased from Sigma- Aldrich. Another cover slip 
is used to cover the top of the inner container with the 
water-decalin interface sandwiched between the top and 
bottom cover slips. Then we overturn the cell. Particles 
fall onto the water-decalin interface caused by the gravity 
for one or two hours and form a monolayer, as shown in 
Fig. 1. The colloidal particles are hydrophilic and will 
prefer to stay in water rather than in decalin. In addition, 
particles keep a certain distance away from the interface 
because of image charge in decalin. The image charge 
force is repulsive for particles approach the decalin- water 
interface from the medium (water) with the higher di- 
electric constant [2J|. The water and decalin layers are 
both 0.8 mm, high. The distance between the cover slip 
and colloidal monolayer is long enough to eliminate the 
influence of cover slips. 

With an inverted microscope Olympus Cool- 
SNAPIX71, the motion of colloidal particles was 
recorded by a digital camera (Prosilica GE1050) of 
14 frames per second. Each image sequence includes 
500 consecutive frames (~ 35.7 sec). Each image is of 
size 176 * 176/i7Ti^. From these image sequences, using 
homemade software, we obtain the particle positions 
and construct trajectories with a spatial resolution of 
60-lOOnm. 



III. EXPERIMENTAL RESULTS AND 
DISCUSSIONS 

Interaction potential U{r) of particles was estimated 
from pair correlated function g{r) of particles, r is the 
particle separation. At dilute limit, g{r) is related to 
U{r) through the function of U{r)/kBT = —ln{g{r)). 
In our experiment, g{r) was calculated by averaging 10^ 
particle positions under area fraction n = 0.04 for PS and 
silica spheres. For even higher area fractions n — 0.08, 
measured g(r) of PS and silica spheres could overlap with 
the formers respectively. That indicates the area fraction 
used could be regarded as the dilute limit. Curves of U{r) 
are shown in Fig. 2. Both particles look like hard spheres 
while r > 2.5d. When r is shorter than 2.5d, an attrac- 
tive potential well appears. For silica particles, the depth 



of attractive potential well is around Q.iksT at r = 1.3d. 
For PS particles, the depth is around OAksT at r = 1.6d. 
What we found accords with the results in Ref. [2^, [2^ : 
like charge particles could attract each other near an in- 
terface. PS particle shows a wider and deeper attractive 
potential than silica, because it has more surface charge 
than silica f26|. 

PS monolayer is located at higher position above the 
oil- water interface than silica because PS sphere has less 
mass density and more surface charge. The average dis- 
tance z from particle monolayer to the interface could 
be estimated by measuring the single diffusion coeffi- 
cient of particles. From the particle trajectories s{t) , 
we calculated the single-particle mean square displace- 
ment (MSD) (As2(r)) = (|s(i -f r) - s(t)|2). The self- 
diffusion coefficient Ds is obtained according to the equa- 
tion (As^(r)) = ADs * T for different particle area frac- 
tion n. Figure 3 shows the dependence of the diffu- 
sion coefficient Ds/Dq on n for PS and silica spheres. 
Dq = ksT /i-nrj^d is the diffusion coefficient of a sin- 
gle particle in water. The solid lines in Fig. 3 are the 
second-order polynomial fitted results. 



Ds/Do = a{l-Pn-^n'). 



(1) 



The parameter a stands for the local viscosity felt by a 
single sphere at dilute limit: the larger a is , the smaller 
the viscosity is. The parameter /3 reflects the strength 
of two body His between the spheres: the smaller (5, the 
smaller two-body His. The fitted values of a, (3 and 7 
are given in Table I. 

The value of a can be used to estimate the distance be- 
tween the interface and the sphere monolayer z. Accord- 
ing to the classical prediction given by Lee et al. [23.l27J. 
the distance from particle's center to the interface z can 
be written as: 



a 16(r/u. + rio){a - 1)' 



(2) 



where a is particle radius, 77^ is the viscosity of water 
and rio = 2.5 cP is the viscosity of decalin at 22.5 °C. 
Substituting the fitted values of a into Eq. 2, we ob- 
tain z = 2.3/z?Ti±0.2/iTO (1.4^m±0.1/Lim) for PS (silica) 
sphere. The calculation is consistent well with the ex- 
perimental observation under the microscope: particles 
keeps 1 — 2/i?Ti distance away from the interface. 

The particles slightly fiuctuate in vertical direction un- 
der ksT. The fluctuation amplitude of silica particles is 
smaller than PS because the former has heavier mass 
density. Within the area fraction region we measured, 
particles near the interface could be viewed as a mono- 
layer. We focus our investigation on the lateral motion of 
particles, parallel to the oil-water interface, in following 
part of the paper. 

Following the trajectory s*(i) of individual particles i, 
particles' cross-correlated motion is obtained via the en- 
semble averaged tensor product of the particle displace- 
ments 12]: 



M,y{r,T) =< A4(t,r)A4(i,T) >,^,-t 



(3) 



Here As^(i, r) = s!^(t + r) — s^(i). i,j are particle in- 
dicates, X and 2/ represent different coordinates, and r 
is the separation between particle i and j (shown as the 
inset of Fig. 4). The off-diagonal elements are uncor- 
related. We focus on the diagonal elements of this ten- 
sor product: Mrr which indicates the correlated motion 
along the line joining the centers of particles (called par- 
allel direction), and Mgg which represents the correlated 
motion perpendicular to this line (called perpendicular 
direction). We found that the measured correlated mo- 
tion Mj-r and Mgg are linear functions of r for small lag 
time r. Thus, the cross-correlated diffusion coefficients 
are defined as D|| = Mrr/^r and D± — Mgg/2T. 

Figure 4 exhibits D|| and D±^ of PS and silica spheres 
as a function of r for different area fraction n. Each curve 
of an arbitrary n, whose deviation was given in the legend 
of Fig. 4, was obtained by averaging at least 10^ particle 
positions. The behaviors of D\\ and D±^ of PS and silica 
monolayer are qualitatively similar. With the increase of 
r, Z)|| and Dx of PS and silica sphere decrease following 
the power law. The specific form of the power law will 
be discussed in the following section. When r > 2.5(i, 
the curves of D\\ (r) and D±{r) for different n almost col- 
lapse onto a single curve respectively, i.e., particle area 
fraction hardly affects the particles' cross-correlated mo- 
tion. When r < 2.5d, however, D±{r) is dependent on 
the area fraction n. With the increase of n, the decay rate 
of D±{r) decreases, i.e., the curves become more flat. 

The mechanism of n independence effect mentioned 
above is in line with Oppenheimer et a/.'s theoretical 
description in Ref. [l^- In the particle monolayer, the 
far- field response of the correlated particle motion mainly 
arises from the momentum diffusion through the 3D sur- 
rounding fluid. Usually the viscosity of 3D fluid suspen- 
sion is the function of volume fraction of particles. In 
our system (j) could be regarded as zero almost, no matter 
how the area fraction n of monolayer changes. The vis- 
cosity of 3D fiuid almost keeps as a constant. Thus the 
far-field 3D His between the particles are independent of 
the area fraction n, so does the cross-correlated diffusion 
D|| (r) and D±{r). Cui et al. [ij] also observed the similar 
concentration independent effects for particles confined 
between two plates, but the reasons to their phenom- 
ena are different. In the confined particle monolayer, the 
fluid momentum is absorbed by the solid boundaries, and 
the far-field fiuid response arises solely from mass-dipole 
perturbation, which is not influenced by the presence of 
neighbor particles. 

The situation will be more complicated for short r 
(< 2.5d), because the particles reveal an attractive poten- 
tial in the range (Fig. 2). Both His and thermodynamic 
interaction involve in particle's correlated motion now. 
For the potential is attractive, the correlation of particle 
motion will become stronger in the range. As a result, 
the curves become more flat. For His, the near-field re- 
sponse of the correlated motion mainly arises from the 
momentum diffusion through the two-dimensional (2D) 
colloidal layer itself. The 2D monolayer viscosity rjm in- 



creases with the area fraction n of the particles. As a 
result, the correlated motion is a function of n. 

To focus on the influences of the oil- water interface on 
the correlated diffusion and eliminate the n effects, we 
average the data of different n in Fig. 4. Figure 5 shows 
the averaged correlated diffusion coefficient I?||(r) and 
Dj_{r) normalized by aDo, which is the single particle 
diffusion coefficient obtained by Eq. 1, as a function of 
particle separation r. At the parallel direction, there is 
D\\{r)/aDo = A^^/r'^^K Here, A\\ is the amplitude coef- 
ficient and A|| = 1 is the decay rate of cross-correlated 
motion. The interface affects only the amplitude A\\ but 
not the decay rate A||. The amplitude A|| of the silica 
spheres is larger than that of PS, which is in accordance 
with the results from single diffusion measurement: the 
fitted value of (3 of silica is greater than that of PS in 
Table I. The phenomenon indicates that the strength of 
His between the particle pair is more strong in the silica 
monolayer than in PS monolayer. The amplitude of A|| 
of the silica spheres is found almost 1.2 times larger than 
that of PS, which number is just equal to the ratio of 
the /3 (0.78/0.62 = 1.2). The same decay rate A|| = 1 for 
both monolayers indicates that the form of His is hardly 
affected by the interface in the parallel direction. 

At the perpendicular direction, D±{r)/aDo ~ 
A±/r^-^. Both the amplitude A± and decay rate Ax 
are functions of the distance z. The decay rates Aj_ of 
the silica spheres is larger than that of PS at short r 
(1.8 vs. 1.5, as labeled in Fig. 4). Different from the 
constant An, the decay rate A^ increases with r. At the 
long r limit, both X± in PS and silica monolayers merge 
into a value of around 2, as shown in Fig. 5. The in- 
set of Fig. 5 shows that the scaled diffusion coefficient 
D±{r) — D±{r)/{a(3DQ) for both PS and silica particles, 
where r — r/z. Both curves are overlapped at long region 
of r > 15, which indicates the effect of the liquid- liquid 
interface tends to be saturated. 

We compare our results with the studies before for 
different boundary conditions: the colloidal spheres dis- 
persed in unbounded 3D bulk, confined by the solid wall 
or at the air- water interface |13l - [la |. We find the cross 
correlated motion of colloidal spheres near oil- water in- 
terface is similar to that of a viscous membrane [1^, ITg : 
A|| = 1 and A^ = 2 at long r limit. For a viscous 
membrane, long r limit means that r ^ Ls, where 
Lg = rjm/irji, is Saffman length 28]. 'qm is the 2D mem- 
brane viscosity and rji, is the 3D bulk viscosity of sur- 
rounding fiuid. For our system, limit of long r means 
that r ^ Max{z,Ls]. Max{z,Ls} is the larger one be- 
tween z and Ls- 

Theory in [l^ shows that the D||(r)/aZ?o does not 
depend on membrane viscosity 7]^. for large r, but 
D±{r)/aDo is dependent on r/m- This prediction accords 
with our results of A|| and A^- Our results suggest that 
the colloidal monolayer near the the interface could be re- 
garded as a 'membrane'. Even though the precise value 
of 2D viscosity ry^ of colloidal monolayer is hard to be 
obtained from our present data, the value of Aj_ (1.8 or 



1.5) indicates the ijm usually is very small. For a 'mem- 
brane' with viscosity rjm, the value of Aj^ will approach 
to 2 at low rim limit and approach to for high rjm limit 
|16| . Precisely, D± (r) will tend to a logarithmic form at 
high rjm limit. Hence Ls should be very small also, and 
the limit of long r usually means r ^ z in our system. 

At first glance, the result that A^ — 1.8 for the silica 
spheres is larger than A_l = 1.5 for PS seems confus- 
ing, which conflicts with the fact that a larger A should 
corresponds to a smaller monolayer viscosity. Since the 
silica spheres are closer to the interface, it is located in 
a more viscous environment comparing with PS. Thus, 
A_L of the silica spheres seems should be less than that of 
PS. The fact is that the viscosity r]{z) felt by the spheres 
is the 3D viscosity of the local surrounding fluid, not the 
2D viscosity rjm of the membrane. The influence of this 
surrounding viscosity r/(z) has been removed by l/(a_Do) 
scaling, where a stands for the effect of local viscosity 
T]{z). The result that A_l of silica is larger than A_l of 
PS stems from His modified by the interface. Oil is more 
viscous than water. The flow fleld induced by the particle 
in water is suppressed by the oil- water interface. His be- 
tween particles decay faster with their separation r when 
the monolayer is closer to the interface, which leads to a 
larger measured A_l for a smaller z. 

Similar to the mechanism of a membrane near a solid 
wall 15], His among the particles in our system trans- 
mit through three paths: I) The 2D flow through the 
monolayer. II) The flow through the fluid layer sand- 
wiched between the monolayer and interface. HI) The 
flow through the upper water bulk or the below oil bulk. 
His through path I is dependent on 2D monolayer viscos- 
ity ijm , while His through path II and HI are dependent 
on 3D viscosity rjo and r]^ of surrounding fluid. 

The cross-correlated motion caused by His through 
Path I is a function of n, because rjm is a function of 
n usually. His through Path I only contributes to corre- 
lated motion when the separation r is within the order 
of Saffman length Ls — rjm/iVo + 'i]w)- In our system rjm 
of the monolayer is very small. With r increasing, the 
weight of His through path I decreases quickly. Mean- 
while, His through Path II contributes to correlated mo- 



tion all the time. So, correlated motion is a function of z 
(see Fig. 5). At long separation limit r » Max{z,Ls}, 
the His through path HI dominate. Since r]o and 77^, 
are hardly changed with area fraction n, the correlated 
motion should be independent of n (see Fig. 4). The 
correlated motion is also independent of z for r ^ z. At 
middle range of r. His through all three paths contribute, 
their relative weights change gradually with r. The de- 
cay rate of the cross-correlated motion shows a crossover 
tendency. 



IV. CONCLUSION 

We investigated experimentally the cross-correlated 
diffusion of colloidal particles near the oil-water inter- 
face. Correlated diffusion coefficient D^^{r) and D±{r) 
are independent of area fraction n unless at the short 
pair separation r. The interface affects correlated diffu- 
sion coefficient D^^{r) and D±{r) distinctly. The interface 
enhances the amplitude of D^^{r), but does not affect its 
decay rate. The results indicate: along the line joining 
the centers of particle pair, the interface modifies the am- 
plitude of the His between the particle pairs, but not the 
form of His. At the direction perpendicular to the line, 
the influence of the interface changes with increasing r. 
The interface enhances the decay rate at short r, while 
such influence tends to be saturated at the long r limit 
(r > 15z, z is the distance between the interface and the 
monolayer). The dependence of the decay rate on r ev- 
idences that the weights of His through different paths 
change gradually. 
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TABLE L The distance 2 from particle's center to the in- 
terface and the fitted values of a, jS and 7 in Fig. 3. Here 
a — Ifim is the radius of particles. 



Sample z/a 



P 



7 



PS 2.3 ± 0.2 
Silica 1.4 ± 0.1 



0.87 ± 0.01 
0.78 ± 0.01 



0.62± 0.12 1.54± 0.61 
0.78± 0.02 0.76± 0.03 



Fig. 5 (color online) The mean correlated diffusion 
coefficient D^^/aDo (solid symbols) and D±/aDQ (open 
symbols) obtained by averaging the data for various val- 
ues of the area fraction n in Fig. 4. The square and 
circle represent silica and PS particles, respectively. The 
inset shows the correlated diffusion coefficient D± (open 
symbols) as a function of scaled inter-particle distance 
r(= r/z) for PS and silica spheres. The solid lines with 
the slope —1 and —2 are the guides for eyes. 



Figure captions: 

Fig. 1. (color online) (a) Schematic view of the sys- 
tem, (b) Optical microscope image of silica spheres {d = 
2.0/im) suspended near the water(up)-decalin(down) in- 
terface at area fraction n = 0.18. 



Fig. 2. (color online) Measured interaction potential 
U{r) for PS (solid circle symbols) and silica spheres (open 
square symbols). Area fraction n = 0.04 for PS and silica 
particles. 



Fig. 3 (color online) Measured self diffusion coefficient 
Ds scaled by Dq as a function of particles' area fraction 
n. Different symbols represent the data for different par- 
ticles. The solid lines show the second-order polynomial 
fitting Ds/Do = a{l — (3n — jn^). 



Fig. 4 (color online) Measured correlated diffusion co- 
efficient I?|| (solid symbols) and D± (open symbols) as a 
function of inter-particle distance r with various values 
of the area fraction n for (a) PS and (b) silica spheres. 
The solid lines with the slope —1, —1.5 and —1.8 are the 
guides for eyes. The inset shows the geometry of mea- 
suring correlated diffusion. 
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